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ABSTRACT 

The author employs a 1 6K linked, FORTRAN II version of the 
problem oriented PACER simulation program written by Prof. Paul 
T. Shannon, Dartmouth College, Hanover, New Hampshire to 
develop a non-linked, FORTRAN IV program which is compatible with 
the GE 625 system installed at the University of Kansas. 

The author then formulates a steam power system model and 
tests several configurations of this system to demonstrate 

(1) the efficacy of the modified PACER program and 

(2) the potential usefulness of such a modified PACER 
program in either graduate or undergraduate chemical 
engineering curriculum. 
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Chapter I 



INTRODUCTION 



General . In early 1965 the National Science Foundation estab- 
lished a project at the University of Michigan to evaluate the current 
and future uses of computers in engineering design, education. The 
final report issued by this project observed that the majority of 
analytical work done in upper-level engineering courses was of a 
routine nature, where the solution procedures are generally brief and 
straightforv/ard and where there is normally only one correct answer. 

For this reason, the report stated, engineering design courses had 
fallen into general disregard and in some instances they have been 
eliminated from the curricula. Rather than take a pessimistic view from 
these circumstances, the report concluded that; 

Such courses can, v/ith the introduction of the computer, 
take on a new luster and make significant contributions 
to the education of our engineering graduates . In the 
past, the engineering student has only had time to 
investigate perhaps one or two possible (and usually 
not very effective) designs. With the new hardware and 
software systems, the student designer will for the first 
time have the opportunity to investigate a large number 
of feasible solutions to his design problem. . . . 

In effect, because of his exposure to a large number of 
designs and his ability to manipulate system configurations 
and design variables easily, the student can gain some 
facets of "engineering experience" which cannot now be 
acquired outside the industrial environment (and which 
would require a much longer elapsed time). Thus, in 
addition to its use as an analysis tool, the time shared 
computer with appropriate design-oriented software 
should enable the student to improve his intuition and 
understanding of the nature of good design: this is the 
raison d'etre of engineering design courses.^ 

One of the most promising approaches to the demand for effective 
engineering design education lies in the recent development of problem- 
oriented approaches to programming languages . The problem-oriented 
approach involves the development of a general processor or executive 
program which is capable of handling a general class of problems v;ith 
a minimal amount of programming effort on the part of the problem 
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designer. This frees the designer to concentrate on the primary- 
objective of system design and analysis. 

Of the problem-oriented simulation programs developed thus far, 
one of the most notable is PACER. PACER was developed by Prof. Paul 
T. Shannon of Dartmouth College and was developed primarily to 
process material and energy balances of typical system design problems 
encountered in chemical engineering. Solutions are processed by a 
modular or unit operations concept with iterative calculations under- 
taken Vv^hen recycle conditions are encountered. 

Objective and Scope of Thesis. This thesis will examine the 
historical development of digital simulation programs culminating with 
the development and some applications of PACER. The essential 
concepts and characteristics of PACER as developed by Prof. Shannon 
will be described and illustrated. The author will further discuss the 
modifications made to PACER to make it compatible with the computation 
system installed at the University of Kansas. Additionally, the author 
will trace the development of a semi-realistic model of a steam power 
plant developed to evaluate the modified PACER program . As a related 
objective, the usefulness of the model in an undergraduate chemical 
engineering education curriculum will be explored. 

Finally, the advantages and limitations of the modified PACER 
program and suggestions for its further improvement and refinement 
will be presented for consideration. 



Chapter II 

HISTORICAL DEVELOPMENT OF DIGITAL 
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SIMULATION PROGRAMS AND 
APPLICATIONS TO ENGINEERING EDUGATION 

Historical Development of Digital Simulation Programs . In itia 1 
attempts toward computer simulation of chemical engineering systems, 
Vy/ere made in the early 1950's and were limited almost solely to analog 
computers. The digital computer at that time required tedious machine- 
language programming and it was slow and expensive for such uses. 

2 

Accordingly, it could not effectively compete with the analog computer. 
(Figure 1 presents a schematic illustration of the historical development 
of simulation programs). 

R. G. Selfridge is generally credited with the development, in 
1955, of the first digital simulation program for use on the IBM 701 

3 

computer . This initial effort was followed during the period from 1956 
to 1958 by the development of a group of digital computers known as 
digital differential analyzers (D.D.A.'s). The D.D.A.'s were designed 
to duplicate the functioning of the analog computer, but proved to be 

4 

exceedingly slow and cumbersome . 

One of the most significant advances in digital simulation 
occurred with the development of procedure- oriented languages, such 
as FORTRAN, in the late 1950's. The availability of these languages 
facilitated the heretofore cumbersome programming task by permitting 
the programmer to write his program in terms closely resembling the 
symbolic representation of problem variables and parameters. A com- 
piler converted the symbolic program into machine language instructions. 
The FORTRAN language additionally was provided with a library of 
functions and subroutines which performed commonly encountered mathe- 
matical operations such as logarithims, exponentation, etc. FORTRAN 
coupled with increasingly spectacular improvements in memory capa- 
bility and speed of operations, where operations are currently specified 
in nano- seconds (billionths of a second), gradually began to assert the 
effectiveness of digital simulation programs. 




Figure 1. 
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The relative advantages and disadvantages of digital and analog 
simulation became more distinct as each system was progressively 
developed. It was readily acknowledged that the analog computer 
enjoyed the advantage of parallel operation which obviated the careful 
attention to logical execution sequence and extensive anticipation of 
alternatives required by the sequential nature of the digital computer. 
Digital computers normally require extensive reprogramming for different 
problem structures. Mathematical integration, a natural operation for 
an analog computer, requires specific and extensive programming on a 
digital computer. On the other hand, the digital computer provides 
certain specific advantages which the analog computer cannot approach: 
greater reliability, memory and logic capabilities as well as a decidedly 
better dynamic range. ^ 

An extensive effort was undertaken to incorporate the best 
features of the analog and digital computer into a digital simulation 
capability referred to as a digital- analog simulator. These simulation 
systems generally consist of block-oriented' subroutines which perfonifi 
or simulate the natural functions of an analog computer such as 
integration and differentation. 

Two representative examples of digital-analog simulators are the 
MIDAS and MIMIC programs. The most significant capability of MIDAS^ 
is its integration function which automatically adjusts the integration 
step to a specified size and combines this feature with a fifth-order 
Milne, predict- correct integration routine. These characteristics 
provide the user with assurance that the integration is being performed 
correctly. ® 

The next effort in digital simulation involved the development of 
equation-oriented, digital-analog simulation programs which did not 
rely upon the block-oriented structure imposed on MIDAS and its 
contemporaries. Chief among this group of programs v^^ere DSL/90^ 
and MIMIC. 

MIMIC, introduced in 1965 by the authors of MIDAS, differs in 
concept from the less flexible block-oriented MIDAS. MIMIC offers 
several other features not comparably available in MIDAS, such as 



a wider variety of available subroutines, FORTMN compatibility and 
intuitive solution capability for implicit algebraic equations. 

A separate class of programming languages has been designed to 
simulate systems characterized by a sequence of discrete events such 
as schedule optimization and resources allocation. Representative 
examples of this type of program are GASP,^^ GPSS,^^ SIMSCRIPT,^^ 
and CSL.^3 

The most significant current advance in the field of digital 
simulation is the development of problem- oriented or application- 
oriented programming languages .which are designed to handle a specific 
class or specific classes of problems. Some of the mors significant of 
this type of program are: 

1. BLODIB^^ and ECAP^^ for signal analysis and circuit 
design in electrical engineering 

2. DYANA ° for vibration analysis in mechanical engineering, 
and 
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3. PACER ' for system simulation in chemical 
engineering. 

PACER, an acronym for Process Assembly Case Evaluator Routine, 
was conceived in 1961 by Prof. Paul T. Shannon while he was in the 
employ of Humble Oil and Refining Company. After extensive modifi- 
cation and development, PACER received its first practical test in 
October 1964 when Prof. Shannon joined with Prof. A. I. Johnson of 
McMaster University, Hamilton, Ontario, Canada and Dr. Norman 
Cooke from Canada Industries Ltd. to undertake a simulation of C-I-L's 
300 ton/day contact sulfuric acid plant in Hamilton. This initial under- 
taking was successfully completed and reported in April 1965. The 
sulfuric acid plant simulation involved the use of 42 equipment, infor- 
mation gathering or parameter setting subroutines which were inter- 
connected by approximately 70 information transmittal streams. The 
model thus conceived involves the solution of approximately 500 
simultaneous equations involving almost 1000 stream variables and 200 
equipment parameters. A complete simulation run required 8 minutes of 
com.puter time on the IBM 7040 computer. 
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In 1965 the staff and the senior chemical engineering students 

at McMaster University joined with C-I-L in the design of a 600 ton/day 

sulfuric acid plant based on metallurgical off-gas. The results of this 

study were presented to C-I-L and representatives of other companies 

in April 1966. Today, more than twelve industrial firms are using 

PACER in the analysis of their processing problems, and thirty-nine 

universities are integrating PACER into their teaching and research 
19 

activities . 

Use of Computers and Simulation Programs in Engineering 
Education . In 1959 the College of Engineering of the University of 
Michigan undertook, under the sponsorship of the Ford Foundation, a 
study of the feasibility of broad scale integration of electronic computer 
use into the educational process. The final report of this study 
recommended that: 

(1) All engineering students be introduced to computing 
techniques early in their college education, preferably 
at the freshman or sophomore level. 

(2) That they develop a proficiency in a procedure- 
oriented language such as FORTRAN, ALGOL, MAD, etc., 
and 

(3) Required upper-level engineering courses incorporate 
the requirement for computer solutions of several engineering 
problems of graded difficulty by the student prior to gradu- 
ation .20 

A disciplinary appendix to the final report entitled "The Use of 
Computers in Chemical Engineering Education" cited the purely 
educational aspects of computers in the education process, stressing 
the salutary effects of precise formulation, definition, and logical 
organization required for computer solution of engineering problems. 

A survey of engineering school deans and department heads 
taken in 1962 as a part of the Project revealed that approximately 30% 
of all engineering students attended schools which had a required 
introductory computer course for undergraduate students. The study 
further indicated 45% of engineering students attended colleges or 
universities which required students in some engineering disciplines 
to attend such an introductory course. A subsequent study conducted 
during the 1965-56 academic year showed that approximately 87% of 
all engineering undergraduate students received some required 
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instmction in digital computer techniques. An analysis of this data 

indicates that in the near future almost all undergraduate engineering 

2 1 

students will attend a required digital computer introductory course. 

Quality, however, apparently has not matched quantity as the 
summary report on "Computers in Engineering Design Education" 
concluded. (See Chapter I). Most of the computer solutions assigned 
in upper level engineering courses are routine in nature. One major 
shortcoming in the field of engineering design has been traditional 
reticence of design instructors to employ sophisticated rriathematical 
analysis techniques in their courses. Further, existing computer 
hardware and programming languages, due primarily to their general 
nature, were not amenable to flexible use in design problems. 

However, two recent developments in the field of computer 
technology hold the promise of a vast improvement in the area of 
engineering design education. The first of these new developments is 
the time-shared computer with remote terminal installations. Those 
facilities will permit a large number of users to have almost instan- 
taneous access to the computer and will permit these users to communi- 
cate with the computer in an almost conversational manner. This time- 
sharing capability, coupled with the development of a new hierarchy of 
problem-oriented computer languages and processors such as PACER, 
permits the visualization of the investigation of a large number of 
feasible configurations of complex engineering problems by a student 

in the course of a semester. In effect, the computer will then take 

22 

on the role of an effective imparter of "engineering experience". 



Chapter III 

CONCEFf , ORGANIZATION AND OPEFVnONS OF PACER 
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Concept of PACER, As previously mentioned PACER was developed 
by Prof. Shannon to provide a highly- flexible, systematic method for 
the simulation and design of complex chemical engineering systems and 
processes. The basic concepts underlying the PACER system may be 
summarized as follows: 

1. PACER is a systematic method which permits a progressive 
development of comprehensive models of chemical engineering systems 
or processes . 

2. The processing or operational sequence and the physical 
configuration of the system or process are included in the input data 
provided by the user and can be readily modified to explore alternate 
system or parametric configurations . 

3. Mathematical models of individual equipments are available 
for use during development of the total system or process. 

4. Earlier work and results may be incorporated into the process 
model with minimal effort on the part of the user. 

5. The successful completion of each analysis contributes to 

the development of a librar\>’ of mathematical equipment models which 

23 

may be used in subsequent analyses. 

PACER conceptually then differs from a procedure- oriented 
language such as FORTRAN. FORTRAN permits a problem solver to 
write solutions for mathematical and engineering problems in a language 
which symbolically approximates commonly- encountered engineering 
and mathematical terms. It does not provide a general or even a specific 
solution method. It's greatest advantage lies in its flexibility — the 
ability to represent many types of problems, ranging from the social 
sciences and law to nuclear physics, with the same language. 

This generality also represents a deficiency of FORTRW. Its 
use, in many cases, requires the complete restructuring of complex 
solution programs for only minor system or problem configurations. 
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Such a requirement can be exceedingly cumbersome to a systems 
designer or analyst v;ho is primarily interested in the results of the 
simulation. 

PACER' was developed to provide to a chemical engineer the 
ability to rapidly and accurately solve varying configurations of complex 
chemical engineering systems. PACER accomplishes this by providing 
a rigidly- structured program which is capable, with minimal effort on 
the part of the user, of solving complex chemical engineering systems. 
No tedious rewriting of the program is necessary for changes in config- 
uration. PACER, then, provides the logical next step in computer 
language design-- macro- procedure-oriented programs. 

It permits a user to systematically describe (by the use of a 
highly- structured data deck) the topology and parametric constraints 
of the system to be solved. Major changes in the actual physical 
configuration of a system can be accommodated with relatively simple 
changes to the data deck. It also provides the user the option of 
specifying or not specifying the order in which equipment v;ill be cal- 
culated. If the order of calculation is not specified, PACER enters into 
a trial and error procedure to determine the solution order. 

The PACER concept of operations may be effectively illustrated 
by the presentation of a simple problem involving (1) a flow mixer, 

(2) a heater and (3) a flow splitter. Figure 2 presents the flow diagram 
for such a system, as well as the major components of the data deck 
which v/ould be used to communicate the system configuration and 
parametric representation to the PACER executive routine. A more 
extensive discussion of the data deck and its components is presented 
subsequently in this chapter and Appendix A. 

Organization of PACER. The major elements of PACER are: 

1. A mainline program which serves as the general processor 
for PACER operations and conducts all iterative calculations. 

2. Information modifying and parameter setting subroutines. 

These include: 

(a) EQCAL - Calls proper equipment siabroutines to calculate 
system or process equipments. 
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a . Flow Diagram 



EQPT. 



STREAM NUMBERS 



1 


MKER 


1 


2 


-3 


0 


0 


2 


HEAl’ER 


3 


5 


-4 


0 


0 


3 


SPLTR 


4 


-5 


-6 


0 


0 



b. Process Matrix 
EQPT. 

NO. EQUIPMENT PARAMETERS 



1 

2 

3 



0 0 0 
0 17. 350. 

0 .2 .8 



0 

2 . 

0 



0 

43.5 

0 



0 

0 

0 



c. Equipment Parameters 
STREAM 

NO. STREAM VARIABLES 



1. 0 15. 

2. 0 15. 

3. 0 15. 

4. 0 15. 

5. 0 15. 

6. 0 15. 



80. 100. 0 
60. 200. 200. 
0 0 0 

0 0 0 

0 0 0 

0 0 0 



100 . 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 



d. Stream Variables 



STREAM VARIABLES ABSOLUTE FRACTIONAL CON VERG ENCE C RI TERIA 
1. 1. 1. .01 .01 .Oi .01 .01 

e. Convergence Criteria 



EQPT. MODE 

1 1 

2 2 

3 3 



f. Equipment Calculation List 
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(b) TEST - if calculation is in an iterative (recycle) mode this 
subroutine tests the values of the output stream variables for the equip- 
ment being calculated to determine if the calculated values meet con- 
vergence criteria established by the user. 

(c) XDATA - provides the space for additional equipment 
variables if required by the user. 

(d) MAJOR - called by PACER when order of calculation of 
equipment is not supplied by user. MAJOR directs search for list of 
equipments which will calculate as output streams up to three assumed 
known input streams. When such a list is developed the equipments are 
calculated as either direct or iterative calculations depending on their 
configuration. 

(e) EQUIP - serves as a bookkeeping subroutine for all equip- 
ments which are calculated directlyo Bookkeeping functions for equipment 
calculated in an iterative mode are performed by the mainline program. 

(f) GUESS 1 - forms under the direction of MAJOR an ordered 
list of unknown input streams. From this list, up to three unknown 
streams are assumed known and a search is made to find a list of equip- 
ment which will calculate as output streams the assumed known input 
streams . 

(g) GUESS2 - carries out under the direction of GUESSl the 
search for a list of equipment which will calculate as output streams 
those unknown streams assumed known in GUESSl. If such a list of 
streams can be found, GUESS2 calls for appropriate calculation. 

(h) DREAD - reads data deck for problem configuration (run). 

A new' set of data is required for each configuration (run) . 

(i) EXREAD - reads additional equipment parameters if pro- 
vided by user as Vv^ell as order of calculation if specified. 

(j) SEFLAG - sets numerical indicators (flags) indicating 
status of equipments and streams provided as data. Flags may indicate 
either that an equipment is used or not used or that a stream is a feed 
stream, product stream, internal stream or not used. 

(k) DPRINT - prints data supplied by user. Primarily used as 
a debugging feature. 



(l) EXPRNT - prints additional equipment parameters supplied 
or calculated by PACER as well as order of calculation specified by user. 

(m) RPRINT - prints results of run conducted by PACER, Results 
are printed at completion of each run. 

3. Equipment subroutines which are mathematical models of 
equipments to be calculated. PACER as supplied by Prof. Shannon con- 
tained forty-two equipment subroutines in the PACER equipment subrou- 
tine library. The user may supply such additional equipment subroutines 
as he may require. 

4. Data deck - the data deck describes the system or process 
configuration and any known stream variables or equipment parameters. 

It also includes certain dimension and control parameters as well as 
certain other operational indicators or controllers. The principal con- 
stituents of the data deck are: 

(a) Process Matrix - the process matrix is the principal 
segment of the data deck. It describes the sytem or problem configuration 
and indicates stream interrelationships. Column 1 contains the equip- 
ment subroutine which performs the mathematical calculation of the equip- 
ment listed in column 1 . The columns after column 2 indicate the numbers 
of the streams flov/ing into and out of equipment number listed in column 
1. All input streams are listed first. All output streams must be preceded 
by a minus sign. 

(b) Equipment Parameters Matrix - contains for each equipm.ent 
listed in the process matrix a list of equipment parameters such as weight, 
volume, power, or heat transfer coefficient. Column 1 indicates the 
equipment number. All columns after column 1 contain equipment para- 
meters for that equipment. The equipment parameters are characteristic 

of that equipment only. 

(c) Stream Variables Matrix - contains stream variables for 
each stream listed in the process matrix. The stream variables represent 
such characteristic items as pressure, temperature, specific heat, total 
flow rate, component flow' rates, etc. Stream variables are consistent 
for each stream. Column 1 contains the stream number. Column 2 is a 
flag indicating the status of the stream. If the flag equals 1, the stream 
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is a feed stream. If it equals 2, it is a product stream. If it is equal 
to 0, it is an internal stream. The remaining columns contain the values 
of the stream variables for that stream. 

(d) Convergence Criteria - contains in each column an abso- 
lute fractional convergence value for each corresponding column in the 
stream matrix. Convergence of an item in the stream variables matrix 
occurs Vv/hen the value calculated for that variable does not vary from the 
current value of that variable stored in the stream variables matrix by 
more than the absolute convergence criterion. The convergence criterion 
for a particular variable location (column) is the same for all streams. 

(e) Calculation Order Listing - if provided or specified by the 
user it indicates the order in which the equipment will be calculated as 
well as the mode (direct or iterative) by which the equipments will be 
calculated. Column 1 indicates the equipment number and column 2 

the mode of calculation. A value of 1 in the calculation mode column 
indicates that the equipment is to be calculated directly. A value of 2 
indicates that the equipment is to be calculated iteratively. If the cal- 
culation mode is 3 the equipment is the last equipment in a list of 
equipment to be calculated in an iterative mode. 

(f) Equipment Subroutine Listing - a listing which establishes 
the correspondence between the subroutine name and subroutine number 
to facilitate calling of proper subroutine for equipment calculations. 

(g) Other control and informational parameters and data such 
as problem dimension parameters, problem control parameters, and title 
of simulation run. These items are specifically defined in Appendix A. 



Chapter IV 

MODIFICATIONS MADE TO PACER 
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Modifications Made to PACER . The PACER source deck supplied 
by Prof. Shannon was impressed on a 2400 foot magnetic tape at 55 6 
bits/inch, 80 characters/record, blocked 1. With the assistance of 
the University of Kansas Computation Center, the magnetic tape was 
converted into an interpreted, punched card deck of approximately 
11,500 cards . 

Detailed analysis of the interpreted deck revealed that: 

(1) The program was organized into eight 16K links as well as 
a starter link and several sample problems. The eight links were 
arranged as follows: 

a. LINK 1 - contained mainline segment for conducting iterative 
calculations, calling and testing subroutines and certain equipment 
subroutines used in the contact sulfuric acid plant simulation. 

b. LINK 3 - consisted of mainline segment which directs search 
for list of equipment when order of calculation is not specified, sub- 
routines for formulation of equipm.ent calculation list, and a subroutine 
for recording values calculated in direct calculations. 

c. LINK 4 - composed of a mainline segment which directed 
input/output operations and subroutines which effected input/output 
operations and stream and equipment flagging. 

d. LINKS 2, 5, 6, 7, 8 - contained remaining equipment end 
information modification and gathering subroutines used in contact 
sulfuric acid plant. 

(2) The linking procedure employed in the program was not 
compatible with the GE 625 system installed at the University of 
Kansas . 

(3) The program was written in FORTRAN II, 

Further analysis indicated that the 11,500 cards could be 
reduced to approximately 2800 by removing all those cards pertaining 
to equipment subroutines used in the contact sulfuric-acid- plant 
simulation as well as those used in the saraple problems at the end of 
the source deck. The cards remaining after this modification constituted 
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the PACER executive routine. However, it still existed in a non- 
compatible linked version and was written in FORTRAN II which was 
not perm.itted on the GE 625 system. 

Accordingly, two additional modifications were required to 
make the PACER executive routine compatible with the GE 625 system. 
First, the program had to be converted to a linked version acceptable 
by the GE 625 and the program had to be rewritten in FORTRAN IV. A 
detailed investigation of the organization of the executive routine and 
the simulation model planned for development indicated that the PACER 
executive and the simulation model could be combined into a non- linked 
program which would not exceed the memory available on the 625 . 

The conversion from a linked to a non- linked version created 
several exacting problems. First, the combination of the linked main- 
line segments into one non- linked mainline program created a statement 
number redundancy. Second, it caused the invalidation of a major 
portion of the linked program transfer logic. This problem required 
extensive effort to resolve. 

The rewriting of the executive routine in FORTRAN TV required 
considerably less effort than the conversion to a non- linked program. 
The major changes required were the rewriting of arithmetic function 
statements and input/output instructions. 
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Chapter V 

DEVELOPMENT OF STEAM SYSTEM MODEL AND 

TESTING OF THE MODIFIED PACER EXECUTIVE ROUTINE 

Obiectives of Testing . The principal objectives of testing 
conducted with the modified PACER executive routine were; 

(1) To verify the validity of the modifications made to the 
original version of PACER; 

(2) To demonstrate the capability of PACER to simply and 
accurately simulate the performance of varying configurations of com- 
plex chemical engineering systems. In this regard it should be 
pointed out that specific care was taken throughout the course of the 
modification of PACER to insure that no alterations were made in 
either PACER logic or method of solution. Finally, 

(3) To provide the experimental and factual basis for evaluating 
the potential usefulness of the modified PACER executive routine in a 
chemical engineering curriculum. This objective was facilitated by 
the steam power plant model developed by the author for testing with 
the modified PACER program. 

Development of Simulation Model. Feasibility or validity testing 
of the modified PACER program presented the contrasting alternatives 
of (1) simulation of a segment of the contact sulfuric acid plant by 
use of the equipment subroutines supplied by Prof. Shannon or (2) 
development of a new simulation model. The second alternative was 
selected since it provided the opportunity for development of a more 
thorough understanding of PACER program logic. 

The simulation model selected was that of a conventional 
steam power system utilizing fuel oil as its energy source. The reasons 
supporting this selection were: 

(1) The familiarity of the author with steam power systems. 

This familiarity was gained by the author from service for two years 
as an instructor in Thermodynamics at the United States Nava] 

Academy, as well as practical acquaintance with steam power system.s 
as the primary means for naval ship propulsion. 
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(2) The general applicability of steam systems to chemical 
engineering . 

(3) The ready adaptability of such a model to the undergraduate 
chemical engineering curriculum at the University of Kansas. 

Accordingly, the author undertook the development of the 
essential equipment subroutines necessary for the approximation of 
a steam power system. As envisioned, the steam power system would 
consist of: 

(1) A Turbine 

(2) A Condenser 

(3) A Cooling Water Circulating Pump 

(4) A Condensate Pump 

(5) A Heater Mixer 

(6) A Feed Pump 

(7) An Economizer 

(8) A Boiler 

(9) A Flow Splitter 

(10) A Superheater 

(11) A Fuel Oil Heater 

(12) A Forced Draft Blower 

(13) A Flow Mixer 

Such a system, it was envisioned, would be amenable to demon- 
strating the ability of PACER to perform iterative material and energy 
balances by the multiple application of individual equipm.ent sub- 
routines . 

The basic system as envisioned and developed by the author is 
presented in Figure 3, Figures 4 and 5 present alternate configurations 
of the basic system. Table 1 provides a stream legend for Figures 3-5. 
Figure 6 presents the process diagrams for the basic system on 
enthalpy — entropy coordinates. Analysis and development of mathe- 
matical models for individual components of the basic system are 
presented in Figures 7-19. 
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Description 

Steam 

Exhaust Steam 
Impurities 
Losses 
Condensate, 

Cooling Water 
Cooling Water 
Cooling Water 
Auxiliary Steam 
Auxiliary Steam 
Condensate 
Auxiliary Steam 
Auxiliary Steam 
Feed 
Feed 

Auxiliary Steam 

Auxiliary Steam 

Feed 

Steam 

Fuel Oil 

Boiler Fuel Oil 

Boiler Fuel Oil 

Boiler Blowdown of Impurities and 
Accompanying Feed 

Fuel Oil 

Superheater Fuel Oil 
Superheater Fuel Oil 
Auxiliary Steam 
Auxiliary Steam 
Air 
Air 

Auxiliary Steam 
Auxiliary Steam 



Table 1 . Steam Power Plant Stream Legend 
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W . - steam flow rate 

steam 

- entering enthalpy 

- internal engine efficiency 

- mechanical efficiency 

hj^g - isentropic expansion enthalpy 
h2 - exhaust enthalpy 



Figure 7 . Turbine 
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Figure 8 . Condenser 
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Figure 9 



Circulating Pump 
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Figure 10. Condensate Pump 
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Figure 11. Heater- Mixer 
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Figure 12. Feed Pump 
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Figure 13. Economizer 
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Figure 14. Boiler 



No mathemalcal calculations take place in 
the splitter. The splitter merely serves as 
a convenient manner for splitting the sat- 
urated ste;\m flow produced in the boiler into 
its various components vjithout needlessly en- 
cum.bering the boiler calculations. 



Figure 15 . Splitter 
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Superheater 
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Figure 17. Fuel Oil Heater 
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Figure 18. Forced Draft Blower 
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Figure 19. Mixer 
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Plan of Testing . The development of a series of approximations to 
achieve the objectives established above required a systematic analysis 
to determine if any of the objectives were complementary or interdepend- 
ent, and, if so, in what manner and in what order. 

An analysis of these relationships revealed that the objective of 
verifying the validity of the modifications made to PACER was the objec- 
tive on which all the other objectives depended. Additionally, it was 
concluded that objectives two and three could be achieved as derivatives 
of objective one. That is, if simulation tests could be developed which 
would establish the validity of the modifications and the accuracy of the 
model, then the efficacy of PACER- and the model and their usefulness in 
chemical engineering education could be determined by derivation from 
those results . 

Accordingly, a series of six simulations were scheduled for testing 
with the modified PACER program.. The first. Run 01, was designed to 
approximate the performance of the basic system with the order of calcu- 
lations specified by the user. The second simulation, designated Run 03, 
involved the same basic system but did not specify the order of calcula- 
tions. The accuracy of the solutions was not of primary concern in these 
tests. Rather, the objective was to ascertain if the modified program were 
processing the data in accordance with the methods described in the docu- 
mentation provided by Prof. Shannon and inferred from an analysis of the 
original PACER. 

Two additional system configurations were tested to demonstrate the 
ability of PACER to provide accurate solutions of varying configurations of 
the same system. The first of the altered configurations, labeled Run 02, 
involved all the components of the basic system except the superheater. 
The second altered configuration consisted of the basic system less the 
economizer. This run was designated Run 04. Both Runs 02 and 04 were 
run with the order of calculation specified and unspecified. The results 
obtained from these runs were compared with hand-calculated values and 
corrections as necessary, and were made in the models until the computed 
values agreed with the hand calcuations for all coses. 

The final results of all runs are presented in Appendix B, v/ith 
extensive diagnostic print-out (a PACER option) provided for Runs 01 and 
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03. A guide for the interpretation of the results of the simulations accom- 
panies the results in Appendix B. 



Chapter VI 

ANALYSIS OF TEST RESULTS 
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Accuracy of Results . Acceptable limits of accuracy vary with the 
nature of the problem being evaluated. Quite clearly, the acceptable 
limits of accuracy in determining the configuration of the critical mass 
of a fissionable material would be much more rigorous than those to be 
imposed on the sytem parameters and configuration of a steam power 
plant. In the latter case, slide rule accuracy is norm.ally acceptable. 
Accordingly, a convergence criterion of .02 (2%) — an approximation of 
slide rule accuracy — was selected for all steam variables other than those 
non-varying items such as stream number and flag value. With this cri- 
terion, all runs converged to solution in no more than five iterations. 

The values of stream variable calculated by PACER were compared 
to hand- calculated values and in all cases were within slide rule accuracy 
of hand- calculated values. 

Effect of Specifying Order of Calculations . PACER undertakes trial 
and error (iterative) solutions in all cases where the order of calculations 
is not specified or where an iterative calculation is specified. Under 
these conditons, the advantage of specifying the order of equipment cal- 
culations was not clear. In Run 02 (without superheating) both the calcu- 
lations with and without the order of calculations specified converged in 
five loops. In Run 04 (without economizer) the calculations with the order 
of calculations specified converged in three iterative loops, whereas the 
calculations without the order of calculations specified converged in four 
loops. Runs 01 and 03 (basic system) converged in five and four loops 
respectively. 

The inconclusive nature of these results would seem to indicate 
that the PACER method of determination of order of calculation is as effec- 
tive in promoting convergence as intuitive designation of the order of 
calculation — or at least no more ineffective. 

Another aspect of this consideration is the comparison of the values 
of stream variables obtained with different orders of calculations. Quite 
clearly convergence values will be expected to vary depending on the 
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logical location of equipment in the order of calculations. While this 
expectation was confirmed, differences were in all cases less than the 
convergence criterion. 

Effects of Changes in System Configuration . Two major changes 
to the basic steam power system were simulated with the modified PACER 
program. The first change involved removal of the superheater from the 
basic system. This change necessitated the use of saturated steam in 
the turbine with the result that for the same turbine power requirement, 
steam flow increased 12%. However, fuel oil flow increased only 1%. 
Condenser surface area increased 6%. 

The second change involved removal of the economizer from the 
steam power system. This change necessitated the addition of the heat 
that would have been added by flue gases in the economizer as the heat 
of combustion of additional fuel oil. This change resulted in a 20% 
increase in fuel oil flow. 

Realism of the Models . The qualitative and quantitative nature of 
the stream variables tended in all cases to verify the realistic nature of 
the models developed by the author. For example, the use of saturated 
steam in lieu of superheated steam would undoubtedly require more steam 
to produce the same turbine power. The greater steam flow would norm- 
ally require more fuel oil. However, the fact that the additional steam 
was saturated and the original superheating load was no longer required 
had an offsetting effect on fuel oil flow. 

Additionally, as was demonstrated in Run 04, the loss of the regen- 
erative heating affect of the economizer would quite clearly require more 
fuel oil for the same power requirement. 

While these results tend to verify the realistic nature of the steam 
pov;er system developed by the author, it should be understood that no 
attempt was made in the development of the model to completely depict 
reality. Rather an attempt was made to provide a semi-realistic model 
which could be refined subsequently. It should be noted in this regard 
that most of the data and mathematical bases for equipmient and sytems 
models are based on reliable data developed from many years in the 
design, development and operation of naval steam power plants. 
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Validity of Modifications Made to PACER. The close conformance 
of the values obtained in the simulations with hand-calculated values, 
as well as the general and specific correspondence of results to format 
and arrangement of data as described in documentation provided by Prof. 
Shannon clearly establishes the validity of modifications to the original 
PACER program „ 



Chapter VII 

AN ASSESSMENT OF PACER 
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Genera l . The author has demonstrated the concept, principles and 
method of application of the PACER simulation program by the develop- 
ment of a complex model of a steam power system, and successfully 
approximating its performance on a modified version of the original PACER. 
The validity of the modifications made to PACER have been conclusively 
established by the accuracy of the results achieved. The realistic nature 
of the model has been tentatively established by the qualitative and 
quantitative nature of the results. 

PACER As A Simulation Device . PACER has been demonstrated to be 
a powerful, flexible and efficient method for simulation of complex con- 
figurations of chemical engineering systems. 

However, certain areas where improvements in PACER method of 
operations or concept could be made were noted and are discussed here 
for edification of potential users. 

(1) The addition or removal of equipment subroutines from PACER 
is cumbersome with present arrangement of calling segments and sub- 
routines and LPHA listing. A more permanent, non-mnenonic name for 
equipment subroutines requiring only a change in subroutine name to con- 
form to permanently established but unused name is one method of allevi- 
ating this problem. This would obviate changes in mainline program which 
become more difficult when executive program is on tape. 

(2) Consideration should be given to establishment of permanent 
STRMI and STRMO matrices to accommodate all streams on a permanent 
basis for a run. This would facilitate multiple calling of stream variable 
and information transmittal in a particular run. However, it could develop 
that the ultimate result would be to trade the present inflexibility of a 
temporary STRMI and STRMO matrix for an inflexible permanent matrix. 

(3) Consideration should be made for incorporating capital invest- 
ment costs into equipment subroutines . This information could take the 
form of a cost which varies with heat transfer area, pressure, temperature, 
etc. Incorporation of this data into the basic structure of an equipment 



subroutine is a relatively simple matter — primarily a function of availa- 
bility of cost information and its relationship to design or operating 
parameters . 

Use of PACER In A Chemical Engineering Curriculum . In evaluating 
the results of the simulation of the C-I-L contact sulfuric acid plant 
using the PACER program. Prof. Shannon remarked: 

It should be emphasized that this study, although it 
involved a unique and very close cooperation with 
industry, was in the main an exercise in the education 
of engineering students. The students were concerned 
with the detailed modelling of specific equipment, yet 
had a very broad exposure to a real industrial plant and 
its over-al] behavior. They learned of the difficulties 
and importance of obtaining good plant data, and of 
writing successful programs of greater complexity than 
previously encountered. The project was admittedly a 
very ambitious one; however the large effort expended 
has placed PACER into the hands of the staff for future 
use in the undergraduate and graduate program. 

The staff generally agree that most of the students have 
benefitted by this exposure to the quite new concept of 27 
information handling and decision making utilizing PACER. 

PACER then has already demonstrated its effectiveness in instilling 
in students some aspects of "engineering experience" which would not 
normally be obtainable outside an industrial environment. This, as 
stated in the University of Michigan Report on the uses of computers in 
engineering education, is the raison d'etre of engineering design courses. 

It is further considered that PACER can be an effective means of 
developing a systematic approach to problem solving by chemical engi- 
neering students. The demands of computer programming in development 
of equipment subroutines and systems concepts requires strict attention 
to problem structuring and formulation. 

Recommmendations . In consideration of the results of the simu- 
lations conducted with the modified PACER program and the conclusions 
drawn above, it is recommended that: 

(1) Consideration be given to the incorporation of the modified 
PACER program into existing undergraduate courses involving basic prob- 
lems in chemical engineering. The steani power system as developed by 
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the author is amenable to immediate use on an individual or term project 
basis, and 

(2) After suitable evaluation of the effectiveness of the initial use 
of PACER that consideration be given to further and more extensive use of 
PACER in graduate and research projects at the University of Kansas. 
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General . The most significant feature of the PACER simulation 
method is its ability to perform realistic and accurate simulations of 
varying configurations of complex chemical engineering systems with a 
minimal effort on the part of the user. The element of PACER which 
makes this capability possible is a highly- structured data deck which 
permits the user to fully characterize system configuration and the re- 
straints or parameters under which it operates. 

Successful use of PACER is 'dependent upon strict adherence to a 
prescribed set of instructions for the formulation, assembly and modifi- 
cation of the data deck. It is the purpose of this appendix to provide 
the necessary instructions in a manner which will permit a user to 
employ PACER' s powerful calculational skills without having to master 
either the logic or structure of PACER. One exception exists to this 
condition. That occurs when it is necessary to develop or modify equip- 
ment subroutines or add them to or delete them from the program. In that 
case a minimal knowledge of the structure and logic of PACER is required 
to insure the subroutines conform to PACER logic and that they are properly 
entered into PACER structure. Both of these aspects will be treated in 
this appendix. 

System Conceptualization . The initial step in preparing a system 
configuration for simulation by PACER is to establish the system concept 
and commit the concept to formalization by preparation of a flow diagram. 
This flow diagram should break equipment assemblies into subassemiblies 
and equipments capable of being mathematically modelled. For example, 
a boiler normally contains an economizer and a superheater as integral 
equipments, and when one normally refers to a boiler it is assumed that 
he includes the economizer and superheater if they are provided. However, 
breaking these elements out and handling them as separate equipments 
facilitates mathematical and system modelling as well as increasing the 
flexibility of the model. Further, the functional components if of suffi- 
cient general nature could be used in multiple applications within a 
single simulation or as a related equipm.ent in subsequent simulations. 
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In developing the system flow diagram all equipments, the per- 
formance of which will be simulated by mathematical models (equipment 
subroutines), should be numbered. No redundancy in equipment number- 
ing is permitted although equipment numbers need not be successive. 
While equipments are being numbered, they should also by named. For 
convenience it would be preferable to assign them a six-character alpha- 
numeric name which is mnemonically related to the equipment function, 
e. g. TURBN5 for turbine and CNDNR5 for condenser. No redundancy is 
permitted in the last 3 characters of a unit name, although individual 
characters within a unit name may be repeated. For example, TUR BN5 
and TUR BN6 are acceptable, while SPL TR5 and SPH TR5 are not. 

After all equipments have been numbered and named, all inter- 
connecting streams should be numbered. No redundancy is permitted 
in stream numbers. However, stream numbers need not be successive 
nor is redundancy between stream and equipment numbers prohibited. 

Equipment Subroutine Format . The equipment subroutine is the one 
facet of PACER that poses a programming burden on the user. If an 
equipment subroutine which satisfies the analytical and calculational 
requirements of an equipment which is to be simulated is not available 
in the PACER subroutine library, the user must develop an equipment 
subroutine which accomplishes these objectives. 

Several challenges confront the user developing an equipment sub- 
routine for use in a PACER simulation. The first challenge is to insure 
that the equipment to be modelled is, in fact, the equipment which should 
be modelled. In this respect a general effort should be made to break 
equipment assemblies down into their functional elements, that is, as 
previously mentioned, a boiler into a boiler, an economizer and a super- 
heater. Mathematical analysis is facilitated and simplified by this 
approach. Further, specific effort should be directed toward incorporating 
generality into subroutines. That is, an economizer is in essence a two- 
fluid heat exchanger. If the equipment subroutine can be written for a 
general two- fluid heat exchanger which is in this application called an 
economizer, then the equipment subroutine can be used in subsequent 
future applications Vv/hero a two- fluid heat exchanger is employed. The 
fact of this generality should be documented into the subroutine by 
appropriate comments . 
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Once the mathematical representation of an equipment has been 
developed it is necessary to translate this into the syntax employed by 
PACER. This task can be mastered quickly by knowledge of a few basic 
concepts of PACER. 

PACER performs equipment calculations by using mathematical 
relationships provided by the user in the equipment subroutine and 
stream and equipment data provided by the user in the stream (SN) and 
equipment (EN) matrices of the data deck (to be discussed further in 
section dealing with the data deck) . PACER calls data out of the stream 
matrix into temporary STRMI (inlet streams) and STRMO (outlet streams) 
for the equipment to be calculated. The inlet and outlet streams for a 
particular equipment are designated by the user in the process m.atrix 
(KPM) . Stream variable information is then ready for use in PACER cal- 
culations. Equipment parameters are called from and returned directly 
to the EN matrix for the equipment to be calculated. Cnee the equipment 
has been calculated the values in the STRMI and STRMC matrices are 
returned to permanent storage in the SN matrix. The STRMI and STRMC 
matrices are then available for temporary storage of inlet and outlet 
stream variables for next equipment to be calculated. 

It is then incumbent upon the user to transcribe his mathematical 
notation into a notation which indicates to PACER where in the STRMI, 
STRMC or EN matrices the variable or parameter to be used is to be lo- 
cated. For example, let us convert the mathematical equation for deter- 
mining turbine stream flow rate into a notation usable by PACER. The 
equation (neglecting dimensional constants) is: 

W = P -J- Wk 

In FCRTRAN this is represented as 

W = P/WK 

However, in PACER notation it would be written as 

STRMI (1,10) = EN (NE,3) / (STRMI(1,5) - STRMC(1,5)) 

(where WK= (3TRMI(1,5) - STRMC(1,5))) 

In the above notation 

STRMI(1,10) - indicates the tenth variable in the first inlet stream 

to the equipment being calculated. It also represents the tenth 
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variable in the stream matrix vector for the stream which is the 
first input stream to the equipment. In the context of this sample, 
this variable represents the stream flow for that stream. 

EN (NE,3) - indicates the third parameter for equipment NE. Equip- 
ment NE is the equipment number for the equipment being calculated 
and will be set elsewhere in PACER. The turbine power rating for 
equipment NE will be stored in this location. 

STRMI (1,5) and STRMO (1,5) - represent the fifth variable in the 
first inlet and outlet streams respectively. In this example these 
variables take on the value of the entering and leaving enthalpies 
of the first inlet and outlet streams to equipment NE. 

The effect of the above calculation is to calculate a value for 
steam, flow for the stream number which is the first inlet stream to equip- 
ment NE. After calculation this calculated value is stored in the SN 
matrix releasing the STRMI and STRMO matrices available for further 
calculations. Since the values of EN(NE,3) and STRMI(1,5) and STRMO 
(1,5) were not modified, the corresponding values in the EN and SN 
matrices will not be modified. 

Figures Al. and A2 . are the equipment subroutines for the mathe- 
matical models developed for the turbine and condenser in Chapter 4 (See 
Figures 7 and 8). It should be pointed out that the assignment of variable 
and parameter locations in the SN, EN, STRMI and STRMO matrices are 
the responsibility of the user. Further, stream variable location must 
be consistent for all streams in the SN, STRMI and STRMO matrices. 

That is location 5 in the SN, STRMI and STRMO matrices will throughout 
the course of a run represent the same variable, but not necessarily the 
same value. For example, in the illustrative example discussed above, 
location 5 represented the enthalpy for all streams but the enthalpy for 
all streams was not necessarily the same value and. some streams did not 
even store a value for enthalpy if it were not relevant to the calculation 
of that equipment. In such a case PACER normally assigns a value of 
zero to that location. 

One additional note of caution is that a common and dimension 
deck identical to that used in the mainline program must be the first item 
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C«»ft**PArER common and dimension dfc< 

DIvENSION STRMl(«,l5),STRMCI(t>,'5), STRMO(«. 1 , STRNC0(5,5 ) , 
lSN(5n,15)iSNC(NO.^)»EN(5n#l5).ENr(?50#55,NECAlL(5 0)»NAMF(30>i 
2NA*‘’E? ( 30 )» NAME3 ( 30 ), NAME 4 ( 30 ), kF. 1 (30 ), KS?M( NO , 3 ), K»IS ( bn )*, KP^( 6 ) , 
3DELS( 15) , DFLSC( *5 ) , <PH( 3o,9) ,KErLHG(30 ) . KSFI AC(50 ) . TI n F , 

4 L pH A (30 > ,KOmPNO( 10 ) , AEN( 4, 40 ) ,NEL IS1 (30 . . KSFPNT ( 16) ,PAPEP( < » 1»3) » 

5^NP&pER(20) I 

CO^'HON STRM!, STRMcI, NIN. STPMO, strhco. nout 

C0^'M0N SNiNSlmAV, SnC, NSTLhX, PELS,- DFlSC.KSEm 

CO^'HCN NEUMAx', FNC, NECLPv * ' 

CO‘-MCN KEFLAG, •''EMAX, Ksn AH, NSMAX. ^'ECALL. NEOUIP 

- •• co^'HOn kph, ^4y^'Ax, n3max, title ■ 

COMMON aEN| LFHA, KOMPNO, MOCCMP# NFL I ST 

COMMON Name , N’AMt:2, NAME 3, NAME4, KRET , KP6T2, KRET3, NOaEN, NOFN 
CO'^HON LIMIT, LIMIT2, LOOP, KRACK, KBaCK2, KVOID 
CO-hON KRUN, kskTS, KCLPA^■, LOOPS, NOGO, KAPDS, KPRINT, KONV 
CO^'MON mchain, mchain, klocki 

- - CO'-MON KOUNTl, K0UNT2, KOHNT3, KOIINT4, KOUNT9 

CO'^MON EN, NE, NLIST, LIST 

CO^’MON KUS,NUSMAX, kps, NPS . - . - 

COMMON MUSI, MUS2, NUS3, KCVl . KCV2* KCV3 

CO^'MON N S liSFD, MF URED, kSprmt, paper, NPAPER, NEXT 
CO*"H0n hFl, NElMAX. WFuELI , UFUEL2, I XRA Y, JXRAY, WBDN 



C»«»«tSU^^ROUT IlE TURBN5 CaLCuLAtFS T^E STEAM FLO’^ REduIRFD TO PRODUCE A GIVEN 
TURbiNF POWER OUTPUT, USE IS MADE OF THE RELATIONSHIP THAT FOR A 
STPaM 1URBINE POWER ^ STEAM FLOW X TUPRINE WORK, WHERE TliRRiNE W'OR.< 

Is MAnUESTEP as the PROnuCT OP tHF change In ENTHALPY OF THE STEAM 
In the TURHINE and the TuRRINF mechanical EFFICENCY. 



no-fnclatupe 

sTRM0(i,5) - outlet steam EN TH AL P Y ( RTU/LB ) 

STRM1(1,p5) - INIET steam F N I m a! PV (BTU/I.R ) 

EN(NE,!>) - Internal 6^’GI^.^ err ic ifncy(PErcEnt) 
EN(NE,6) ^ mechanical EFF I C I EMC Y (PEpCrNT ) 

EN(NE,4) - isfnthopic steam EyHAuST ENTHAL P Y ( b TU/LB ) 



♦ •••#cai.cula |E OutUFt steam Enthalpy 

. STRH0(i,5)cST«MI( 1 ,5 )« (STRMI ( 1 . 5 ) - FN ( NE , 4 ) ) • ( E N ( ME . 5 ) / 1 0 0 , ) 
CULATt; STEAM FLOW ... - - 

ST■HMl(l| 8 )c(EN(NF, 3 )* 25450 n 00 ,)/^EN(NE. 6 )•EN(NE, 5 )^^(STRNI(l, 5 )- 

1 EN(NF, 4 ))) - 

STRMl(l,lO)sSTRMI(l,8) 

STPmO( 1 , 8 )bSTRmI ( l.d) .... 

STshO(l,iO)=STRMO(lia) 

STRHOn ,9) = 0, 

STwmI(Ji9)=sTHM0(1*9) 

DO 10 Jtl 1 , 15 . - 

STRhOd, J) 30 , 

ST3MK1, J)aSTRMO(l, J) - 

?.0 CO^TINl'£ 

DO 20 1.6,7 — • - 

STRmO( 1, I )cQ, 

ST®MI (1 , naSTRMOd, I ) ~ 

20 CO^TINue 

C*<^*^.IF MCHAIN • 3 ^ EOUIPMEMT Is BEING CALCULATED DIRECTLY 
IF (HCHAin- 3 ) lOOOilOOliloOO 

1000 Ixrayco . ..... 

GO TO 1002 

1001 IXPAY 33 

1002 RETURN 

END • 



Figure Al, Subroutine TURBN5 
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SUBROUTINE C^JDNR 5 



C 

C* 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

cr 

c 



••••PAfER common and dimension DECK 

DIMENSION STPMI ( B# 15 ) , STRMC I <5 » 5 > , STRMO( 8 , 15 ) » STPMC0(5 » 5) , “ — ’ '* 

l5Nr'50> ! 5 ) # SNC (50 I 5 ) , Esj< 30 » 1 5 ) » ENC ( 30 » 5 ) I NECALL < 50 ) . name < 30 ) » 

- 2NA-^E2t30) »NAMb3(30)#NANE4(30) »KE1 (30) ,K5EM(5n,3) ,KlJS(5n) »KP$<6K ” 

30ElStl5) >DELSC(5) .KPN <30. 9) , KER.AG ( 30 ) . KSEL AG ( 50 > , T I TIE ( , 

4LP^-A(30 ) i;tOMPNO(10 ) . A En < 4 , 4 0 ) . NE L I ST ( 30 , < ) . KSPRNT ( I 6 ) , P ACER < 4 , 15 ) . ^ *" 

5NPApER<;?fl) 

CO^'MON StPHl. STRHCI. NtN.‘" STRMO ; 'STPHCO . NOUT 

C0“MCN SN,NSLMAX, SNC, NSCLMX. PELS. nELSC,KSEM 

CO‘"MON NELMAX, ENC.’NECLHX " 

CO^^HON KErUAG, ^'E^tAV, KSFLAG, NrMAX, NEC ALL > NEOU I P 

CO'’MON KPN, NPMAX, N3MAy, TITLE ■' 

CO'*MCN aEN, L^HAi KOMPNO, NOCOMP, NFLIST 

CO”MCN NAME. NAME2, NAME3. NAME4, K RE T , ' K RE T 2 T K RE T 3 N 0 AFN . NOEN 

COMMON LINIT, limit?, LOQR, KBACK, KBACK2. KVOIP 

CQ-’MON KRUN, KSETS, KCLEAM, loops, NOGO, KAPDS,~KPRINT7“KGNV : 

CO"*MON NCMaIN, ncmaINi KLOCKI 

■ - CO-'MON KOUNTl, K0UNT2)- K0UNT3, K0UNT4 KOimS 

COMMON FN, NE. NLIST. lIST 

CO^'KCN KUS,^’USMAX, KPS, NPS 

common NUSl, MUS2, NUS3. KCVl. KCV2. KCV3 

C0'*MCN NSUSEO, NEUSEP, KSPRMT , PAPPR. NPAPER, next -- --- 

COMMON KEl.NElMAX.WFU'^Ll, t*TUEL2. IXRAY.JXPAY.WMDN 

****SU‘*RCUT INE CNDNR5 DEtER'MNES TM£ COOLING WAtER FLO-J RATE AMD THE sjp- 
FA''E area required to comdemre the exhaust steam from. the turbine, 

IT also calculates the flow PATE OF IMPURITIES INTROnUCEH INTO THE 
SYRTE.x and THE STEAM AMD WATER LOSSES RY LEAKAGE, E \/A P OP A T T DM , . G T C . 

th- condenser surface area is cai.culated fpqm The relationship that 

C r U X A X THETAM, VHERC li 1$ TMF OVERALL HFa^ TRANSFER COEFriClcVT 
and THFtam is THE LOO MEAN TFMPEpaTURE DIFFErEnCF. 

NO“EMCLAtuRE 

STRmO( 1,8) - Lossrsf LR/HP) ' 

STRMI(1,6) - TURBINE FXHAUST R TF A M ( 1.0 /HR ) 

Factor for DTERMI^Mr^lO lorses(lb lossrs/lr steam) 

!N3UPlTlF3(Lf^/HR) 

EN(.AE,3) - factor for DfttRVI^'IMO I MPUR I T I ES ( lR I NPUR I T I ES/LS STEAM) 
STRMI(.;,3) - ^AKE-U° FEFD(L0/nR), EOIJAL^ LOSSES ♦ BLOWDOWN PLOW RATE 
COOLING water flow ratE(L8/HP) --- -- - 



EN ( Nc , 4 ) 
STrmi ( 2 , 8 ) 



STR^ I (3, b > 



C* 



C* 



- — iFacop,Eo.n GO to 77 — 

Go TO 78 

•77 UBRNrO. - ■ • 

7b Co^TINUE 

»ca'.cllate losses ' ■ 

ST-N0(l,ft)=STRMI<l,8)*FN(NF.4) 
ST-MO(1,10 )sSTRmO(1,B) - ■ - 

«^T»*CAi culate Impurities 



■ s T^M I ( 2 , 8 ) =Er; ( ME , 3 ) *STRm I < 1 , 8 ) '' ’ ‘ 

ST-M0(2,11)=STPMI (2,8) 

STPMI (1,10)=STRMI (1,8) 

ST .-Ml (2, ii)=STPM! (2,8) 

c#****cai.culate make-up feed ■ ' 

STRMI (4, 8)=STRM0(1»8 )+W90n 

STNMI (4,9)=STRNI (4,8) * * • 

c#***#cai*cllate heat transferred 

CsfSTRMI(l,8)-STPM0(l,8))*(STRMla, 5 ) - S TRMO ( 2 , 5 ) ) ♦ STRM I ( 2 . 8 ) # 
' KS.TRMI (2.5)"^TRM0(2,5) )^STRMI (4,8)« (STRMI t 4 , 5 ) - S TRMQ ( 2 , 5 ) ) 

C#e**f calculate cooling water flow — 

STFMl (3,8) =Q/ (EN(NE,5 )« (StPM0(3. 4)-StRM I ( 3. 4) ) ) 
STPM0(3,8)=STRr:I (3,8) ' 

STPMl (3, 15)=STPM I (3,0) 

: '• STPMC(3,15)aSTRND(3,8) 



Figure A2. Subrouilno CMi3MR 



calculate log mean TEMf'E^ATURE DlFFEREr'CE 
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in each subroutine. This facilitates communication of information from 
subroutines to mainline and other subroutines. 

Incorporation or Removal of Equipment Subroutines into PACER 
Structure . The flexible nature of PACER permits the incorporation or 
removal of equipment subroutines to conform to core memory limitations, 
to permit the simulation of completely different systems or to permit 
modification of existing systems. Such flexibility is not achieved, how- 
ever, without a commensurate demand on the user. In this case, it is 
the requirement that the incorporation or removal of a subroutine requires 
a change in the semi-permanent structure of PACER. 

Each subroutine, in addition to being assigned a six- character, 
alpha-numeric title, is assigned an equipment subroutine number depending 
on the location of its calling statement in PACER. Equipment subroutines 
1-10 are called from EQCAL and subroutines 11-30 are called from the 
mainline program. The location of these calling statements and maximum 
number of subroutines are a function of the system to be sim.ulated. How- 
ever, a note of caution. An increase or decrease in the maxim.um number 
of equipment subroutines or any other dimension establishing parameter 
requires a change in the appropriate matrix or vector in the common and 
dimension deck. 

The incorporation of a subroutine requires 

(1) Obtaining an equipment subroutine number from those not 
currently being used. 

(2) Writing a calling statement for that subroutine. 

(3) Providing appropriate transfer directions or transfer indicators. 

For example, the addition of an equipment subroutine COOLR8 as 
subroutine number 5 would require the insertion of the following FORTRAN 
cards in EQCAL immediately following statement number 4 and its transfer 
card, if any, and immediately preceding statement number 6: 

Column No . 5 7 

5 CALL COOLR8 
RETURN 

The insertion of the same subroutine as equipment subroutine number 15 
would require the insertion of the following cards in the mainline program. 
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immediately following statement number 14 and its following transfer 
card: 



Column N o. 4_ 7 

15 CALL COOLR8 
GO TO 7777 

Or the addition of COOLR8 as equipment subroutine number 21 or any 
number greater than 21 would require the insertion of the following cards 
in the proper location following statement number 5201; 

Column Noo 4 7 

21 CALL COO LR8 
GO TO 7777 

The removal of an equipment subroutine requires the reverse procedure — 
removal of the calling statement and the transfer card and replacement of 
these cards with the ti'ansfer card assigned the statement number of the 
deleted subroutine. For example, if equipment subroutine 15 were deleted 
the FORTRAN cards described above would be replaced by the following 
FORTRAN statement: 

Column No. 4 7 

15 GO TO 7777 

One additional change is required in addition to those decribed 
above when an equipment subroutine is incorporated into or removed from 
PACER. This change is necessary due to a feature of PACER which estab- 
lishes the correlation between equipment subroutine names and numbers 
and permits the multiple calling of equipment subroutines for different 
equipments. That is equipments number 5 and 10 may both be pumps 
which can be calculated by subroutine PUMP5. This feature referred to as 
the LPHA listing stores the last three characters of the equipment sub- 
routine name in the position corresponding to its equipment subroutine 
number. The format of the LPFIA listing, 12/ (15A3, 10A3), permits the 
inclusion of the first 25 equipments on card i. If 30 equipments were 
used two cards would be required — 25 on the first card and 5 in the first 
15 spaces of the second card. If an equipment subroutine number is not 
occupied, its three characters are left blank. A control card indicating 
the number of equipments for v/hich an LPHA designation should be read 
precedes the LPHA listing and is indicated in the format statement as an 
12 format. 
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If equipment subroutine C00LR8 v;ere assigned subroutine number 
15, the following entry would be made in colum.n 43-45 of the first card 
of the LPHA listing: 

Column No. 43 - 45 
L R 8 

If equipment subroutine PUMP5 were designated equipment number 28 the 
following entry would be made in columns 10-12 of the second card of the 
LPHA listing: 

Column No. 10 - 12 
P 5 

If, on the other hand, equipment subroutine PUMP5 were removed from the 
PACER structure, columns 10-12 of the second card would be left blank. 

Data Deck . The data deck is the means by which the user conveys 
to the PACER executive routine is a detailed discription of the system to 
be simulated. Only one simulation run may be processed for each data 
deck supplied. However, several simulation runs may be processed 
consecutively by providing consecutive data decks for each run. 

In general, a data deck for a simulation run is composed of 15 
groups of data. However, the first group, the LPHA listing, is provided 
only with the first data deck of a multi-deck execution. All decks after 
the first deck consist of only 14 data groups. The data groups and their 
required order of assembly are listed and discussed below. In addition, 
data decks for the six simulations run to test the feasibility of the modified 
PACER program are provided as illustrative examples at the end of this 
appendix. (See Figures A3. - A8.). 

LPHA Listing . Format: 12/(1 5A3 , 10A3) . This data group indicates 
the allowable names of PACER equipment subroutines and provides the 
requisite correspondence between the equipment subroutine number, the 
equipment subroutine name and the equipment number. The first card 
contains in 12 format the num.ber of equipment names to be read from the 
LPHA listing. The remaining cards contain the LPHA listing of the last 
three characters of the allowable equipment subroutines. No redundancy 
is permitted in the three- character names although individual character 
redundancy is permitted. The absolute location of an equipment subroutine 
in the LPHA listing must coincide with the number assigned to the equipment 
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subroutine. That is, the equipment subroutine name placed in the 25th 
location in the LPHA listing m.ust be the last three characters of equipment 
number 25 . 



Title Card . Format: 2(12A3). A single card provided for the assign- 
ment of a unique identifying title to each run. 

Control Parameters . Format: 1015. A list of 8 parameters which 
control the processing of PACER simulations. The parameters, in their 
order of listing on the input record, are: 



1. KRUN - Run number. 

2. KSETS - Controls intermediate printing of PACER status. 
Following is a listing of allowable values and options each 
provides . 



Value Option 

0 Prints initial conditions of matrices , inter- 

mediate conditions of matrices every KPRINT 
loops and final results of run. 

2 In addition to above, prints values in STRMI 

and STRMO after each calculation and num.ber 
of stream variable deviations from convergence 
criteria . 



3 In addition to above, prints flag lists for 

streams and equipments. 

3. KCLEAN - If set equal to zero, it causes PACER to zero 
all matrices and vectors before reading in data. 

4. LOOPS - Maximum number of loops to be permitted in 
iterative calculations. 

5. NOGO - Controls simulation options if convergence is 
not obtained in number of LOOPS pennitted. If NOGO = 0, 
the answers obtained during the last loop will be taken as 
correct and the calculation will continue. If NOGO 7^ 0, 
additional trial and error testing of assumed input streams 
will be made until either convergence is obtained, execution 
time expires or calculation is determined to be impossible. 

6. KARDS - If greater than zero, provides a deck of cards 
punched with the stream variables matrix. 

7 . KPRINT - Controls number of loops between printout of 
stream variables matrix. 

8. KONV - Not used. Leave blank. 



Dimension Parameters . Fonrnat: 1015. This data group contains 

nine parameters which control the allowable dimensions of a run. Any 

change in these variables requires an identical change to the parameter 

where it appears in the dimension deck. The following is a listing of the 

nine parameters in the order of their recording on the input record: 

1. N2MAX - maximum number of rows pennitted in the process 
matrix. 
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2. N3MAX - maximum number of columns permitted in process 
matrix. 

3 . NEMAX - maximum number of equipments permitted in 
system to be simulated. 

4. NSMAX - maximum number of streams permitted in system 
to be simulated. 

5 . NELMAX - maximum number of columns permitted in EN 
matrix. 

6. NECLMX - maximum number of columns permitted in equip- 
ment control matrix. 

7. NSLMAX - maximum number of columns permitted in SN 
matrix. 

8. NSCLMX - maximum number of columns permitted in stream 
variables control matrix. 

9. NOCOMP - maximum number of components in each stream. 
Not used other than for information purposes. 

Process Matrix (KPM ) . Format: 15, 1X2A3, 1X2A3, 1X915. The 
process matrix contains the topology or configuration of the system to be 
simulated. It is preceded by a single card (Format 15) which indicates 
the number of rows to be read in the process matrix. Each rov>^ corresponds 
to an equipment in the system configuration. Equipment numbers need not 
be consecutive but must be in ascending numerical order. The first item 
to be read from the input record (Format 15) is the equipment number. The 
second and third items are the first throe and second three characters, 
respectively, of the 6 character equipment subroutine used to calculate 
that equipment. The next 9 items are the input and output streams for that 
equipment. Input streams are listed first and output streams, preceded by 
a negative sign, follow. A word of caution is necessary at this point. One 
must remember that equipment subroutines are written in anticipation of or 
in conjunction with the order of input and output streams. Accordingly, 
any changes in the order of input and output streams for an equipment must 
be checked against the equipment subroutine to insure that revised input/ 
output stream listing is consistent with equipment subroutine. If a dis- 
parity exists either the equipment subroutine or the order of input/output 
streams must be revised to bring them into conformity with each other. 

Equipment Parameter Matrix (EN) . Format: 5FI5.5. The EN matrix 
is preceded by a single control card (Format 15) which indicates the number 
of equipment parameters vectors to be read in the EN matrix. The first item 
read in an equipment parameters vector is the equipment number. The 
second item on the input record is not used and is left blank. The remaining 
thirteen items are the equipment parameters for the equipment in question. 
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The equipment parameters are unique for each equipment and normally 
represent such parameters as power, heat transfer coefficient, etc. The 
equipment parameters for each equipment employed in the steam power 
system are listed in Figure A9 . 

Equipment Control Parameters Matrix (ENC ) . Format: 5F15.5. The 
ENG matrix is preceded by a single card (Format 15) which indicates the 
number of rows in the ENC matrix. Representative data to be used in the 
ENC matrix would normally include allowable differentials across equip- 
ments, costs, etc. Control vectors are only included for those equipments 
which employ them. If no control vectors are employed, the control card 
should be punched zero or left blank and input/output should proceed to 
next data segment. 

Stream Variables Matrix (SN) . Format: 5FI5.5. The SN matrix is 
preceded by a single card which indicates number of stream vectors to be 
read into program. If control variable is zero or blank, the input/output 
will proceed to next data group. The stream variable matrix contains such 
variables as pressure, temperature, enthalpy, total steam flow, etc. 

This location labeling is consistent for each subroutine. The SN vector 
for the steam power system is structured as follows. 



STRM. 


FLAG 


PRESSURE 

(PSIA) 


TEMPERATURE 

(°F) 


ENTHALPY 

(BTU/LB) 













SPEC. WT. 
(LB/FT^) 


VELOC. 

(FPS) 


TOTAL FLOW 
(LB/HR) 


H2O FLOW 
(LB/HR) 


STEAM FLOW 
(LB/HR) 













INPURITIES FLOW 
(LB/HR) 


AIR FLOW 
(LB/HR) 


FUEL FLOW 
(LB/HR) 


FLUE GAS FLOW 
(LB/HR) 


COOMNT FLOW 
(LB/HR) 
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Stream Control Matrix (SNC) . Format: 5FI5.5. A single card 
precedes SNC and indicates number of stream control variable vectors to 
be read. If this value is zero, PACER proceeds to next data group in data 
deck. The SNC matrix contains control parameters such as viscosity, 
volume, etc. The distinction between the variables specified for this 
matrix and that provided for SN are arbitrary and lack discrimination. 
Accordingly, SNC need not be used unless required. It is not employed 
in the steam power cycle. 

Stream Variables Test Vector (DELS ). Format: 5FI5.5. This group 
is preceded by a single card (Format 15) which indicates test vectors is 
to be read from DELS matrix. If control card = 0, no test vector is to be 
read and PACER skips to next data group. If control card value ^ 0 then 
PACER will read one test vector NSLMAJC long. This same test vector 
will be used as test vector for determining convergence of each stream 
variable for all streams . 

Stream Control Variables Test Vector (DELSC) . Format: 5FI5.5. 
Preceded by single control card which indicates if test vector is employed. 
If control card value = 0 (Format 15) then no control test vector is employed 
and PACER will proceed to next data group. If control card value 7 ^ 0, 
then test vector is employed and is NSCLMX long. Only one test vector 
is employed for stream control variables for all streams employing control 
variables. This test vector will check convergence of control variables 
for all streams employing control variables. 

Preferred Stream Numbers (KPS ) . Format: 1015. KPS is preceded by 
a control card (Format 15) which indicates number of preferred streams. 

If control card value equals zero, then PACER proceeds to next data group. 
If control card 7 ^ 0, PACER reads numbers of preferred streams, i. e. those 
streams which will be assumed known first in a trial and error calculation 
to determine order of calculations. Normally these should be recycle or 
other internal streams since feed and product streams are assumed known 
in trial and error calculations . 

Extra Equipment Vectors (AEN) . Format: 5FI5.5. AEN is preceded 
by a control cord. If control card value (Format 15) equals zero, PACER 
will proceed to next data group. If control value 7 ^ 0, PACER reads into 
AEN extra equipment parameters for number of equipments indicated on 
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control card. This matrix is only used if insufficient space is available 
in EN and ENC. Length of vector is arbitrarily established at 40. 

Order of Equipment Calculation (NELIST) . Format: 1015. NELIST 
is preceded by a control card which indicates number of equipments con- 
tained in NELIST matrix. If the order of calculation is not specified, 
the control card value (Format 1015) equals zero and PACER proceeds to 
next data group. Each card consists of five pairs of data. The first 
value of a pair represents the equipment number. The second value 
represents the calculational mode for that equipment. If mode indicator 
equals one, the equipment will be calculated directly. If the mode 
indicator equals two, the equipment will be calculated as an equipment 
in an iterative list of equipment. Finally, if mode indicator equals three, 
the equipment will be calculated as last item in list of equipment to be 
calculated in an iterative mode. 

Selected Streams Vector (KSPRNT ). Format: 1615. Contains as the 
first item the number of stream variable vectors to be printed each itera- 
tion. If it is not desired to print any stream variable vectors each itera- 
tion the card should be left blank. The stream numbers for which the 
stream variables are to be printed are listed after the first item. A maxi- 
mum of 15 stream variable vectors may be printed each iteration. 
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Format and Interpretation of Results . As discussed in Chapter V, 
six simulation tests were developed to test the validity of the revisions 
made to the original PACER program, the efficacy of PACER as a simulation 
device and the potential usefulness of PACER as an instructional device 
in a chemical engineering curriculum. The results of those simulation 
tests are presented as Figures Bl. - B6. 

The final values of the stream variables matrix and the equipment 
parameters matrix are presented for each run. In addition for runs 01 and 
03 extensive diagnostics of the PACER solution process are presented 
for analysis by the user. This is a PACER option. The stream variables 
are consistent for all streams in the stream variables matrix. That is, 
stream variable 5 will represent enthalpy for all streams. However, the 
enthalpy for all streams will not necessarily be the same. 

For the steam power system developed for this project, the following 
stream variable arrangement was decided on and should be used as a 
guide in interpreting values of stream variables in the stream variables 
matrix. Note that there are five columns in each stream. This dimension 
may be parametically varied by the user by adjusting dimensions of SN 
matrix. 
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10.0000 

0 , 

0 . 

11 .0000 

62.4000 

б, 4969 

12.0000 

0 , 

0 , _ 

13.0000 

0, 

0 , 



14,0000 



0 . 


500,0000 


600.0000 


1344 .9000 


0 , 
0, 


169936,5604 

0. 


-- 0. 

0* 


16993 P .5604 ■ 
0, 


0, 


2,0000 


126.0000 


1011. 2000 


0, 


169930,5664 — 


o ; 


169930.5604 “■ ■■ 


c . 


0, 


0 . 


0. * . ’ 



1.0000 15,0000 80.0000 49.oogo 

O', 8,4 969 0. o; 



0 , 




0, 


0, 


2.0000 




15,nono 


126 .booo 


94 ,0000 


0 , 


509,8157 


0. 


509.8157 - 


0 , 


0 . 


0, 


0 . 


0. 


2, OOno 


126 .0000 


94 . 0000 


0 .- 


171646,449? 


171637.9512 


0. ■ ■ 


0 . 


0 , 


0 . 


0 . 


0 . 


30 ,nonn 


70. 0000 


0 . 


0 . 


6403132, noon 


__ . 


0 , ■ 


0 . 


n, 


0. 


6403132.0000 


2.0000 


30 ,0000 


95.0000 


0. 


0. 


•* 6493132, noon 


0. 


0. 



0. 0, 0. 6403132.0000 



1.0000 

0. 

0. 


15,0000 

6403132,0000 

a. 


70 .0000 
0. 

0, 


0. 

0 . 

6403132.0000 


0. 
0 i 

0, 


5CO.00O0 
" 2041,6575 

0. 


467.0000 

0. 


12o5.00'10 
2041.6575 

0. 


0, 

0. 

0. 


30,0000 

2041,6575 

0. 


250.0000 

0. 

0. 


1058 . 7000 
-7041.6575 

0. 


0, 

0. 

0. 


30 .0000 

171646,449? 

0, 


126.0000 
1 71637.9512 

0. 


94.0831 
0 . ’ 

0 . 


0, 


500,0000 


467 . 0000 


1205.0000 


0, 

0. 


104 ,7822 ■ 

0, 


0. 


10^ . 7822 

ft. 


0, 

0* 

0. 


30,0000 

104,7822 

0, 


250.0000 

0. 

0. " 


1 058.7000 
104.7622 
0. 


0. 


30,0000 


162.0036 


130.0867 
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6? .40no 


0. 


176569, ?773 


1 76560 . 779.3 


0. 


e ,4969 


0. 


ft. 


0. 


n . 


15 . nono 


0. 


5on , nono 


162.0036 


• 131.4815 


6?.40f*0 


0. 


176569 , ?773 


1 785^0 . 7793 


ft . 


e, 4969 


0. 


0. 


0. 


ft . 


16 ,fl000 


0. 


50ft .noon 


467 . nnoft 


1?05 . ftOftO 


p» 


0. 


1829,781? 


0. 


1829.781? 


p . 


0. 


ft. 


0. 


0. 


17 . none 


ft. 


3n , nopft 


25ft. POOP 


1 n58 , 7pno 




0. 


1829 . 781 2 


0 . 


182*?. 7812 


p . 


0 . 


ft . 


0 . 


ft . 


IP .oono 


0. 


500 , ftpon 


367. noon 


3 .3 9 , ft 0 ft 


n . 


0* 


176569,2775 


178560.7793 


ft , 


6.4969 


0, 


0 . 


ft. 


ft . 


1 ’ . 0 0 0 0, 


ft. 


5 ft 0 , ft ft ft 0 


467. OftOft 


1 ?ft5 . ftOftft 


Oc 


0. 


173858,4316 


ft . 


173858.4316 


n. 


0. 


0 . 


0 . 


ft , 


2P . pono 


0. 


3 ft . n 0 ft 0 


140 .OOOft 


ft . 


P « 


0. 


9238,1 2?8 


ft . 


ft . 


p. 


0. 


9238 , i22« 


0 . 


ft , 


21 ,6000 


1 .OftOn 


1 5 , n r f> n 


7 ft . 0 0 ft ft 


ft , 


c. 


0. 


91^8, 4 nP8 


P . 


ft . 


0. 


0 


9l^P,4ftP8 


0. 


0 . 


2 ?,pono 


0, 


•3 ft , 0 0 n ft 


140,0000 


ft . 


ft. 


ft. 


9288, 1228 


0 . 


ft . 


ft. 


0. 


9788 , 1 2?8 


ft . 


ft , 


23, none 


2 .onon 


5 0 ft , ft ft ft 0 


4 6 7 , ft 0 P n 


449 , nnno 


ft. • 


0. 


1707.8826 


16‘’9.38S7 


ft . 


6,4969 


0. 


ft. 


fl. 


ft . 


2<.0000 


0. 


3 ft , ft P ft ft 


1 4 0 , 0 n 0 0 


ft . 


n . 


0. 


1455.3858 


0 . 


ft . 


0. 


0. 


1455 . 38^8 


0 . 


ft . 


25,0000 


0 . 


3 P , ft 0 ft n 


1 4 0 . ft p 0 ft 


ft . 


0. 


0. 


1455 , .3838 


0. 


ft . 


ft . 


0. 


1455,3838 


0. 


0 . 


26,0000 


i.onon 


15 , ftfton 


70.000ft 


ft . 


ft. 


0. 


1455, .3838 


0. 


ft. 


0. 


0. 


1455, .38^8 


«0 . 


n . 


27,0000 


0. 


/ * 

3ft , nnm 


335. noon 


1 705 . onne 


0, 


' 0. 


361 .346,3 


0 . 


361 . 3460 


0. 


0. 


ft . 


0 . 


0. 


26,0000 


0. 


3ft . ftftft ft 


?5p .nnoft 


21 9. nnno 


c. 


0. 


361 . 


0. 


361.3460 


0. 


0. 


ft . 


0. 


0 . 
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25,oono 

0 , 

o» 

30.0000 
0, 

0. 

31.0000 

0 . 

0 , 

3? , 0000 
0. 

0 , 

33.0000 

0. 

0. 

3^ * 0000 

0 , 

0. 

35.0000, 

0 , 

0| 

36 . 0000 . 
0 . 

0. 



1,0000 

0 , 

17340<1.2363 



0 I 

300.0000 

1734P4.2363 

0. 

0 , 

0. 

0. 

0. 

0. 

2, ooon 

0 , 

0 . 

0 , 

0. 

0. 

0. 

0. 

0 , 

1 . 0000 

0* 

0 , 



15 , nooo 

173484.2363 
0^ _ 

15.0000 

173484.2363 

0 , 

500 .0000 
2616,0356 

0 . 

30 .0000 
261^,0356 

0 , 

15 . 0000 

If 4227.7422 

0. 

500.0000 
169938.5604 

0 , 

15,0000 

184227, 7422 

0 , 

15,0000 

2207,2014 

0 , 



80 .0000 

0. 

80, .0000; 

0,. 

Or . 

467.0000 
0 . 

0. 

' 250,0000 

0. 

0 . 

2006,4997 

0 . 

1 84227.7422 

467 . 0000 
0 . “ 
0. 

3000,0000 

0 . 

184227.7422 

70,0000 
2209 . 2014 
0. 



0 . 

0 , 

jn_._ 

0 . 

0 . 

0 , 

1205.0000 . 

0 . 

1079,6000 
2616 , 0356 

0 , 

0 . 

0. ■ 

n, 

1205.0000 
169938 . 5684 . 

0. 

0 . 

0, ' 

0 , 

38 . 0000 

0 . 

Q . 



CONDEi^StR surface A‘>EA REOUIRFD i 



5650, 75S0, Ft, 



ECO^iOHI7ER SURFACE AREA REQUIREO 



55.58SQ. FT. 
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STR£Ah variables 
0. 

— 0, 

0. 


matrix at run 

500,0000 

16993S,56B4 

0, 


— 


, 


1,0000 


■'3 'LOOP 

630.0000 
■ “ ■ 0. 

0. 


1344.0000 

169933.5634 ^ ^ 

0 , 


0 1 
0| 


2,0C00 


0, 


2,0000 


1P6.0000 


1011,2000 * 


* 0. 


0. 


169933,5684 


0. 


169938,5634 


0 • 


0. 


0. 


0. 


JJ! J.. 



3,0000 1,0000 15,0000 BO. 0000 49,0000 






8,4969 


0 . 


0. 


0 . 


0 . 


4,0000 


2.0000 


15.0000 


126.0000 


94,0000 


■ — 0 , 


0 , 


509.8157 


0 . 


' 509,8157 


0 , 


0 , 


0 , 


0 . 


0 . 


5.0000 


0. 


2,OOno 


196.0000 


94.0000 


_ 62^4qo0 


0 . 


■ 171646,4492 


171637.9512 


0. ----- 


6,4969 


0. 


0 , 


0 . 


0. 


6,0000 


0 , 


30.0000 


70.0000 


0. 


0 , ' * 


0 , ' ~ 


‘6403132,0000 


0 , * 


0 . 


^ 0 , 


0 . 


0 , 


0. 


6403132,0000 


7,0000 


2.0000 


30,0000 


95.0000 


0. 


- -• 0 . 


0 . 


6403132,0000 


0 . 


■ — 0 , 


0 . 


0 . 


0 , 


0. 


6403132,0000 


6,0000 


1.0000 


15,0000 


70 .0000 


0 . 


'■64,0000 


0. 


6403132,0000 


. 0 . 


0. 


0 « 


0 , 


0, 


0. 


6403132.0000 



9,0000 


0 


, 500.0000 


467.0000 


1205.0000 


jj ^ . 


0, 


r 2041,6575 


or 


“2041,6575 


0, 


.9.' 


, 0. 


0. 


0. 


10,0000 


0 


, 30,0000 


250,0000 


1058,7000 


0 ■ 


0. 


, 2041,6575 


0 , 


2041.6575 


0 . ^ 


0 


. 0 , 


0, 


0. 


11,0000 • 


0. 


, 33.00QO 


126.0000 


94,0831 


62,4000 


0, 


171646,4492 


171637.9512 


fl. 


6,4969 


0. 


• 0, 


0^ 


0. 


12,0000 


0. 


, 500,0000 


467,0000 


1205,0000 


0 , 


0 . 


, 104,7822 


0 , 


- 104 .7622 ■ 


0 , 


0. 


. 0 , 


0. 


0. 


13,0000 


0. 


. 30,0000 


250.0000 


1053,7000 


0 , ' 


0 , 


, ■ 104,7822 


0. 


“ 104,7022 


0, 


O' 


. 0, 


0 . 


0, 


1 14.0000 


0. 


, 30,0000 


162.3491 


130.4322 
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62.4000 0, 176639,2773 176630,7793 0, 

6,4969 0 , 0 , ' 0, “ ■ 0. ~ 



5 , 0 000 C . loo ,-OOno~" 16273491 131^6270 

62,4000 0 , 176659,2773 178630,7793 0 . 



: 6,4969 - 

16,0000 ' - 

. 0, 


0, 

0, ^ *” 

Of 


■ 0, 

500 ,0000 

1830,4984 


0. 

467.0000 

0. 


120570000 

1830.4984 


0 » 


0 


0 , 


0. 


0 , ” 


17,0000 


0 , " ' 


30,0000 


250 ,'OOOC 


1058.7000 - , 


0. 


0, 


1830;4964 


0. 


1830,4984 


0, - 


0, 




0. 


-- 07 '■ 


16,0000 


0 . 


500,0000 


"367 70000 


33970000 


0. 


0. 


176639,2773 


i7«63o; 7793 


0. 



8,4 969 - ” 0, ^ 0, 0. 0. 



19.0000 0, ““ 50 0,0 00 0 4 67,0 000 1205.0 000 



0. 


0. 


176861 , 3945 


0. 


176861.3945 


— 0," 

- - - 20,0000 


0, 

0, 


0, 

' 3o,oono 


’ 0. 

140.0000" 


-- 0, 


0, 


0. 


9448,3510 


0. 


0, 


0, 


0, ■ ~ 


9448,3510 


0, 


- 


21 ,0000 


1,0000 


15,0000 


; 000 n “ 


‘ 0 , 


0. 


■ 0, 


9288,1228 


0. 


0. 


;o, - 


0. 


“ 9288,1228 


0, 


0, ■ ■■ 


— 22,0000 


0, 


30 , 0000 


— 140.0000 " 


^ 


;_0i 


0. 


9448,3510 


0. 


0 . 


0, 


.0, 


9448,3510 


0. 


° * 


23,0000 


2,0000 


500,0000 


— 467,0000 


449,0000 


0, 


0. 


1707,8826 


1699,3857 


0 . 


- - - 8,4969 

< 


0. 


0 , 


0, 


' " 0 , ■ 


24 ,0000 




30,0000 


_ 140.0000' 


**0, ■ ■ — 


0, 


0. 


1455,3838 


0. 


0. 


0, 


0. 


~ 1455,3838 


0, 




25 ,0000 


0 , 


“ 30,0000 


■ 140,0000 


0 • "" 


0, 


0. 


1455,3838 


0. 


0 , 


' ® » 


0 , 


1455 , 3830' 


0, ~ 


Q_._ 


■ 26,0000 


1,0000 


‘15,0000 


70,0000 





0, 


o« . 


1455,3036 


0, 


0. 


0, - “ 


0, 


1455,3038 


- 0. 


■ “• 


27,0000 " ■ ' 


0 , 


30,0000 


335,0000 


1205, nooo 


1 

0. 


0. 

0, 


366,1073 
p, 


0. 


366,1073 

■ 0. - 


28.0000 


0. 


‘ ‘^O.OOOO 


■ ‘ ■ 250,0000 


219,0000 


0, 


0. 


366,1073 


0. 


366,1073 


0 , * 


0 , ' 


0 , 


0. 


0, 
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29,0000 


l.ooon 


15. noon 


80 , OOOO 


0, 


0, 


0, 


175770.1719 


0. 


0. 


0, 


173770.1719 


0, 


0. 


0. 


3o,nooo 


0, 


15,0000 


80 , OOOO 


0. 


0. 


300. OOnn 


175770,1719 


0. 


0. 


0, 


175770.1719 


0, 


0 . 


0. 


31 . OOOO 


0, 


500 , OOOO 


467.0000 


l2o5 . OOOO 


0. * 


0. 


2650,5062 


' c. 


“ 2650,5062 


0, 


0. 


0 , 


0 . 


0 . 


3? .0000 


0. 


30.0000 


250 . OOOO 


1079,6000 


0. 


0. 


2650,5067 


0. 


2650.5062 


0, . 


0. 


0, 


0. 


0. 


33,0000 


2,0000 


15,0000 


2020 .7875 


0 . 


0, 


0. 


106673,9062 


0. 


0. 


0 , 


0, 


ft. 


1 86673.906? 


0 . 


3^,0000 


0, 


500 , OOOO 


467 , OOOO 


1205,0000 


0. 


0. 


169938,5664 


0, 


- 169938.5684 


0, 


0. 


0. 


0. 


0, 


35,0000 


0. 


15,0000 


3000.0000 


ft. 


0, 


Ot 


186673,9062 


0. 


0 , 


0. 


0, 


0. 


1 86673.9062 


0 . 


36,0000 


1.0000 


15,0000 


70 . OOOO 


38 . OOftO 


0, 


0. 


2209,7014 


2209.2014 


ft. 


0. 


0. 


ft. 


. 0. 


0. 



CONV^PGEKCt= OBTAINED ATTFR LOOP NUM0f:ff A 
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Ei^UIPHENT parameters MATRIX 



- PINAL 



1,0000 
90,0000 
0 , ' 



0 , 20000,0000 926,0000 60,0000 

0 , ‘ ‘ ~ 0 . 0 . 0 .^ 

C * 0 , 0 , 0 , 



2 ,0000 

650,0000 

0 , 

3.0000 

0 . 

c, 

4.0000 

0 , 

0 , 

5.0000 

0 , 

0 , . 

6.0000 

0 , 

0 , 

7.0000 

n , 

n , 

8.0000 

0 , 

0 , 

9,0000 

0 , 

0 . 



0 , 

0 . 

0 , 

0 , 

0 , 

0 . 

0 . 

0 , 

0 , 

0 . 

0 . 

0 . 

0 . 

c . 

0 . 

0 . 

0 , 

0 , 

’ 0 . 

200.0000 

0 . 

0 * 

0 . 

0 . 



0,0001 


0,0030 


0 ,9700 


0 , 


0 , - 


0 . 


_ 0 _. 


0 . 


0 . 


70,0000 


93 .0000 


996,0000 


0 , ■ 


■ 0 , 


0 , 


0 , 


o .< 


0 , 


70 ,0000 

ft — 


93.0000 


996,0000 


^ f 


.0 . 
0 . 


0 , 
0 , 


1 ,0000 


0 . 


0 . 


0 , 


0 . 


0 , 


0 , 


0 . 


0 . 


70 ,0000 


93,0000 


996 .0000 


0 , 


0 . 


0 . 


0 , 


^ 0 . 


0 , 


30 O .0000 


1,0000 


0.2000 


0 , 


0 . 


0 . 


0 . 


0 . 


0 . 


0,8800 


0 . 1200 


0 . 


80 ,0000 


449 ,4000 


0 . 


0 . 


0 . 


n . 


0 , 


0 . 


0 . 


0 , 


0 , 


0 , 


0 , 


0 , 


0 , 



10 ,0000 

0 . 

0 . 

11.0000 

0 , 

0 , 

12,nooo 

49.0000 

996,0000 

13.0000 

0 , 

0 , 



0 , 

0 , 

0 , 

0 . 

0 . 

0 , 

0 . 

32, 0000 

0 , 

0 , 

0 , 

0 , 



0 , 6800 

80,0000 

0 , 

0,4800 

0 , 

0 . 

60,0000 

26 ,0000 

0 . 



0 . ' 
0 . 



. 0,1200 

0 , ‘ 

0 . 

0 . 

^ 0 , 

0 , 

95.0000 

115,0000 

- 

0 , 

0 , 

0 , 



c. 

0 . 

0 . 

0 , 

0 , 

0 , 

4. 0000 

109,0000 

0 , 

0 , 

0 . 

8 . 
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STREAM variables HATRIX ' RUN 








"3 - HNAL 




1.0000 


0. 


900.0000 


680,0000 


1344,0000 


0 . 


0 t 


169938 ,5684 


0 . 


169938 , 5684 


0. 


0. 


0. 


0. 


0. 


2,0000 

... 


0. 


2.0000 


126,0000 


1011.2000 




0 » 


169938 ,3684 


0. 


” 169938,5684 


0. 


0. 


0, 


0 * 


8 . 


3,0000 


1.0000 


19,0000 


80,0000 


49,0000^ 



- n, 0 . 8,4 969 or 0. * " 

8.4969 0 . 0 , 0 . 0 . 



4,0000 


2,0000 


15,0000 


126,0000 


94,0000 


“ 0 ; 


0 I 


'509,8157 


0. 


509.8157 


0. 


0, 


0, 


0. 


0. 


5,0000 


0. 


2,0000 


126,0000 


94 , OOOO 


^ 62.4000 


„ 0 . 


' 171646,4492 


^171637,9512 


0 ^ 


6,4959 


0. 


0. 


0. 


0, 


6,0000 


0. 


30 , 0000 


70.0000 


0 . 


— 0, 


” * 0. 


6403132 ,0000 


0. 


0 , 


0. 


0, 


0 . 


0. 


6403132 ,0000 


7,0000 


2.00C0 


30 ,0000 


95,0000 


0 . 




O. "^ ~ 


' 6403132,0000 


0. 


0 , 


: 0 , 


0, 


0. 


0. 


__ 6403132.0000 


8,0000 


1,0000 


15,0000 


70,0000 


0. 


" 64,0000 


0. 


6403132,0000 


0. 


- 0, 


0, 


0. 


0 ^ 


o._ 


6403132.0000 


9,0000 


0, 


500.0000 


467,0000 


1205,0000 


0 ( 


0 . 


2041,6575 


or 


204176575 


0 . 


0 . 


0. 


_ 0^ 


0 . 


10,0000 


0. 


30,0000 


250,0000 


1058 , 7000 


- 0 . 




2041,6575 


0 ^ 


2041.6575 


. 0. 


^ 0. 


a. 


0. 


0 , 


11.0000 


0, 


30 , nono 


126.0000 


94,0531 


62.4000 


0, 


-171646,4492 


171637.9512 


0. -- 


8,4969 


0. 


0. 


0. 


0, 


• 12,0000 


o; 


500,0000 


467,0000 


1205, OOOC 


Q ^ - 


— 0. 


104,7822 


0. 


104.7822 


0. 


0. 


0. 


0. 


0, 


' 13,0000 


0. 


30,0000 


250.0000 


1058.7000 


- 0 , 


0 . 


10 4,7822 


0 • 


I04 , 7822 ~ 


0. 


0. 


0, 


0. 


* 0. 


14,0000 


0. 


30 , oono 


162,3491 


130.4322 


62,4000 


0, ■ ' 


“178639,2773 


178630.7793 


■■ ■ — ' 0 , 


8,4969 


0. 


0. 


0. 


0. 


15,0000 


0. 


900.0000 


162.3491 


131.8270 
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6?. 4000 


0. 


178639,2773 


178630.7793 


gi 


6,4969 


0. ' ■ 


0 • 


o; 




16.0000 


0, 


500,0000~ 


467,0000' 


12C3,GOnO ■ 


0. 


0. 


1830,4904 


0. 


1830,4984 


«• 


, - 0. 


- “0. ' 






17,0000 


- 0, 


30,0000 


250,0000 


1058.7000 


0. 


0. 


1830,4984 


0. 


1830,4984 


0 1 


0. 


0,' '■ 


0. 




1«,0000 


0. 


500,0000 


367 ,0000 


339,0000 


0 • 


Ot 


178639,2773 


176630,7793 


0* f ‘ 


8,4969 


- -- 0." 


' ■ ■ 0, - 


0. 


0, 


19,0000 


0, 


500,0000 


467 .0000' 


1205,0000 


0,. 


0. 


176861,3945 


0. 


176861 .3945 


■ 0, 


0. 


“ 0. 


■ 0, ■ ‘ 


0, 


20,0000 


• 0, 


30,0000 


140,0000 


" 0^ -- 


0. 


0. 


9448,3510 


0, 


0, 


' 0 • ■ “ 

21,0000 


0 1 

- 1,0000 


9448,3510 

15,0000 
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